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We study how dispersions of colloidal particles in a cholesteric liquid crystal behave under a 
time-dependent electric field. By controlling the amplitude and shape of the applied field wave, 
we show that the system can be reproducibly driven out of equilibrium through different kinetic 
pathways and navigated through a glassy-like free energy landscape encompassing many competing 
metastable equilibria. Such states range from simple Saturn rings to complex structures featuring 
amorphous defect networks, or stacks of disclination loops. A non-equilibrium electric field can also 
trigger the alignment of particles into columnar arrays, through defect-mediated force impulses, or 
their repositioning within a plane. Our results are promising in terms of providing new avenues 
towards controlled patterning and self-assembly of soft colloid-liquid crystal composite materials. 


PACS numbers: 61.30.Jf 81.05.Xj 81.16.Dn 47.57.jd 

Liquid crystal colloids are attracting considerable at¬ 
tentions for their potential as novel and versatile ma¬ 
terials. These include photonic crystals, self-quenched 
glasses, meta-materials, new biosensors and multistable 
devices m- At the same time these composite materi¬ 
als are of considerable fundamental interest. In fact, col¬ 
loidal objects disrupt the orientational order of the host¬ 
ing liquid crystal generating a remarkable array of topo¬ 
logical excitations, including simple Saturn rings (§}, ne¬ 
matic braids mm, and exotic knotted or linked discli¬ 
nation networks [Mig. While these structures may 
form with spherical particles in nematics or cholesterics, 
complex defect patterns can also be achieved by dispers¬ 
ing colloidal particles with nontrivial topologies such as 
handle-bodies, ribbons or knots [16l - [l9] . 

Intriguingly, these defect textures can, in turn, medi¬ 
ate interactions between colloidal particles [9I 120H221 . and 
create an inherently glassy free energy landscape. The 
associated elastic energy barriers typically dwarf thermal 
noise and hence trap the system in a variety of different 
metastable states. This property is key to applications, 
because it opens the possibility to build novel devices 
capable of switching between different soft metastable 
states, each having specific optical or mechanical proper¬ 
ties. Because the switchable states are long-lived, these 
devices are multistable and exhibit memory effects, which 
are desirable features for the design of new-generation 
energy-saving technologies such as e-paper (23fj27| . Ac¬ 
cordingly, devising novel colloid-liquid crystal composites 
with defect properties that can be tuned by external in¬ 
tervention, represents a major standing challenge of soft 
matter physics. 

Towards this goal, here we explore a new way to create 
robust distinct defect and particle patterns by exploit¬ 
ing a time-dependent electric field acting on a cholesteric 
hosting electrically-neutral colloids. As shown by the 
studies in refs. [281129] . time-dependent electric fields hold 
great potential for being used as powerful manipulation 


tools alongside with well-established techniques based on 
constant fields or optical forces EB EB. The behaviour 
resulting from the competition between the nonequilib¬ 
rium driving mechanism and the free energy landscape 
complexity is very rich: even for a single colloid, the 
time-dependent field can create a variety of defect struc¬ 
tures which depend sensitively on the details of the ex¬ 
ternal force. When several particles are dispersed in the 
cholesteric, this external tunability can be harnessed to 
drive a reproducible global rearrangement of the colloids 
from random to columnar or planar arrays. The results 
can have practical implications for the realization of com¬ 
posites with controllable optical, rheological or mechan¬ 
ical properties. 

The system that we consider consists of a set of spher¬ 
ical colloids embedded in a cholesteric liquid crystal, and 
subject to the action of a spatially-uniform, but time- 
dependent, electric field, E. Accordingly, the Landau-de 
Gennes free energy density, /, results from the linear 
superposition of the bulk, elastic-distortion and electric 
field terms, / = fb + fei + Je- Adopting the repeated 
indices summation convention, their expressions are I32j : 

fb = t (! - 5) Qlfi ~ *FQa,fiQfrQ<ra + (Q 2 afi ) 2 , 

fel — ~2 3) 3” (Sfl H - czfi') ] ? 

Je 127r ^a-^aQ '/3 (1) 

where Q is the traceless and symmetric tensor order pa¬ 
rameter [321 S3], A 0 sets the scale of the bulk free energy 
density, 7 controls the magnitude of the ordering, L is 
the elastic constant, p = 2ir/qo is the cholesteric pitch, 
e a is the (positive) dielectric anisotropy of the material 
and e a /3 7 is the Levi-Civita antisymmetric tensor with a , 

/? and 7 indexing the Cartesian components. 

To match the properties of typical cholesterics, which 
have pitch ss lpm, free energy density scale ~ 10 5 Pa and 
elastic constant ~ 6 pN we set the units of time and 
length to lps and ~ 0.03/im respectively, and fix A 0 = 1, 
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7 = 3 , L = 0.065 and q 0 = 2n/32. The correspond¬ 
ing dimensionless chirality k = ~ 0.3 

and reduced temperature r = 27(1 — 7 /3)/7 = 0 [32] 
are such that the cholesteric phase is the lowest free en¬ 
ergy phase at zero field [311 EE] • The degree of coupling 
with the latter is captured by the adimensional parame¬ 
ter 8 = ^27e a E2/(32irA 0 'y). Here we choose the val¬ 
ues 8 \ = 0.0377, 8 2 = 2\/28^ , £3 = 48±, that span 
a range wide enough to straddle the bulk cholesteric- 
nematic transition, see Supplementary Material (SM) for 
a numerical estimate of the onset of the transition [36]. 

The phenomenological equation of motion for Q is [58] 

Qa/s = ~ |Tr(fq)<W) , (2) 

where T = f v f is the total free energy and the collective 
isotropic rotational diffusion coefficient, T is set equal 
to 0.5, corresponding to a rotational viscosity of ~ 1 
Poise. Eq. [2] is solved numerically with a finite difference 
scheme on a grid of size N x = N y = N z = 96 with 
periodic boundary conditions. The off-lattice positions of 
the centers of the spherical colloids (with radius Ro = 8 
corresponding to 0.25 /mi), x, evolve according to an 
overdamped Langevin dynamics, 

*« = F a /£ + y/(2 k B T)/£Ut), (3) 

where £ a is a Gaussian white noise of unit variance, the 
Stokes friction coefficient is £ = 6niqR 0 and the effective 
viscosity, 77 , is taken equal to the rotational one. F is 
the force acting on the particle which comprises a steric 
and a liquid crystal mediated contribution (see SM [36] 
and [58]). The homeotropic anchoring at the colloid 
surface is enforced by setting Q a p = q[n a np — 6 a ,p/3) 
where n is the unit vector normal to the surface and q 
is proportional to the magnitude of order (q = 1/3 is 
appropriate for 7 = 3 ). As a first, reference case we 
discuss how switching on and off an electric field affects 
the texture of a cholesteric hosting a single colloid. The 
system was first relaxed in zero field in its most stable 
state where the cholesteric helix axis is along the z di¬ 
rection and the particle is embraced by a -1/2 twisted 
disclination line m E2i- Next, an electric field pulse 
along z was applied from time t = 0 to t = 0.24s, see 
Fig-0 The pulse duration is long enough to allow the 
system to relax at all three field strengths 81 , 82 and £3, 
see SM [36]. For each case, Fig. [I] illustrates the com¬ 
mon initial condition, the held-on relaxed configuration 
and the post-pulse asymptotic state. At the smallest 
field, 8 \ , the field-relaxed configuration is only slightly 
distorted (with a modest 20% contour-length increase) 
with respect to both the initial and post-pulsed configu¬ 
rations, which are identical. However, both the interme¬ 
diate field 82 and the largest one £3 are strong enough 
to induce a nematic ordering of the medium so that the 
twisted disclination straightens up into a Saturn ring de- 



FIG. 1: Evolution of an equilibrated cholesteric colloid (top) 
subject to an electric field pulse starting at t = 0 and ending 
at t = 0.24s. The asymptotic states at t = 0.24s for the three 
increasing field strengths, £1 < £2 < £3 are shown in panels 
Ai, Bi and Ci, respectively. The corresponding post-pulse 
relaxed states at t=0.42s are shown in panels A 2 , B 2 and C 2 . 


feet in the field-relaxed configuration. Apart from this, 
the system behaves differently at the two fields. First, 
the nematic order induced by £3 is practically perfect, 
while it is only partial for £ 2 , where several extended 
disclination lines perpendicular to the field direction are 
observed. Second, these disclinations reflect the presence 
of several metastable states which render the relaxation 
dynamics at £ 2 substantially slower than either at £1 or 
£ 3 , (see SM [56]'). 

These differences become even more pronounced after 
the field pulse, when the liquid crystal texture cannot re¬ 
lax back to the equilibrium cholesteric order because it is 
trapped into metastable structures that are created dur¬ 
ing the post-pulse relaxation and whose qualitative fea¬ 
tures depend on the pulse strength (i.e. the system has 
memory of its history). In particular, after the strongest 
field (£ 3 ) is switched off, the Saturn ring surrounding 
the colloid expands, while two planar hoops normal to z 
appear at the colloid surface close to the poles. At the 
same time, a regular array of stacked disclination loops 
appears above and below the colloid at a vertical spacing 
about p/2. These defects appear with the same rela¬ 
tive colloid positioning in larger simulation cells, unlike 
the small disclinations on the vertical edges in panel C2 
which, being finite-size effects, move further away from 
the colloid, see SM. Compared to £ 3 , the relaxation dy¬ 
namics for £ 2 follows a qualitatively-different pathway. 
While the defect evolution in proximity of the particle is 
similar to the £3 case, the network of defects away from 
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it is much more complex and arrests in an amorphous 
structure. The post-pulse asymptotic free energies of the 
states subjected to the fields £2 and £3 exceed the initial 
one by respectively, ~ 2.5-10 i ksT and « 2-10 4 &bT. Such 
large offsets are indicative of glassy-like landscapes where 
a plethora of energy barriers much larger than thermal 
energy prevents relaxation to the lowest free-energy state. 

As we discuss below, the rich morphology of the free 
energy landscape, which is both created and exposed by 
the transient action of the external field, can be exploited 
to stabilize several metastable states characterized by dif¬ 
ferent symmetries and different responses to an external 
perturbation. This can pave the way for the realization 
of switchable devices based on cholesteric colloids, and 
controlled by time-dependent external electric fields. 

To illustrate this possibility, we considered an electric 
field modulated as a square wave with strength £3 and 
switching time r = 0 . 02 s, that is longer then system re¬ 
laxation time, see SM [55] . 

Fig-0 illustrates the time behavior of the system free 
energy of the system over a period of the square-wave 
field. The superimposition of the free energy relaxation 
curves for different field switching events reveals an al¬ 
most perfect periodic steady-state response of the sys¬ 
tem which hops between the two asymptotic field-on and 
field-off states shown in Fig. [l] in panels Ci and C 2 , see 
SM [55]. This periodic response is also observed in the 
dynamics of the defect structures whose typical config¬ 
urations within a single oscillation period are shown in 
Fig.[2]panels A 1 -A 4 . Note that throughout the cycle, the 
colloid position is practically fixed and the Saturn ring 
expands remaining planar. 

Next we considered the steady-state system response 
to a sinusoidal wave. Again, the colloid is essentially 
motionless during the cycle and is surrounded by a sin¬ 
gle Saturn ring at the peak field strength, when the 
medium consistently attains a purely-uniaxial state. Re¬ 
markably, however, the sine-wave modulation leads to a 
significantly different evolution of defect patterns as the 
field intensity decreases. In fact, the colloid surface be¬ 
comes transiently encompassed by three hoops. Further¬ 
more, the periodic pattern of loops which bridge various 
cholesteric domains away from the colloid now grows into 
an intricate defect network before dissolving to leave a 
nematic state when the field returns to its peak strength 
(see Fig.J 2 J.B 4 ). The above phenomenology demonstrates 
that a time-dependent field provides a promising avenue 
to reproducibly navigate a system of cholesteric colloids 
through its complex landscape of metastable states. The 
remarkable underlying physical behaviour and its prac¬ 
tical ramifications are aptly exposed by extending con¬ 
siderations from a single colloid to a collection of par¬ 
ticles. In this case, cooperative effects are expected to 
arise from effective inter-colloid interactions mediated by 
the intricately-evolving defect patterns la [ 22 ]. 

Accordingly, we considered a dimer oriented along the 
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FIG. 2: Steady-state evolution of a single cholesteric colloid 
in a periodically-modulated electric field. The cyclic time- 
evolution of the free energy and the field are shown in the first 
and third panel from the top for a square-wave and sine-wave 
modulation, respectively. For each case, four representative 
snapshots of the system dynamics are shown. 


x direction and relaxed in zero field. The equilibrium 
state features two unlinked and unknotted disclination 
lines which wrap around the particles (Fig. [ 3 ] and mm 
thereby producing an effective colloid-colloid attraction 
of ~ —4 • 10 2 ksT. A square-wave electric field, along z, 
with period 0.02s and strength £3 was next switched on. 

The typical subsequent evolution of the system is il¬ 
lustrated in Fig. [3] When the field is switched on for the 
first time the cholesteric undergoes a transition to the 
nematic state, associated with the formation of a figure- 
of-eight entangled state around the dimer mm- When 
the field is switched off, the system remains trapped in a 
metastable state, reminiscent of the one found for a single 
inclusion (Fig. [TJ panel C 2 ). After this, however, the in¬ 
teraction between the two colloids gives rise to a different 
and unexpected evolution. In particular, after two cycles, 
the colloids displace, mostly along z and, when the field 
is switched on again, the figure-of-eight defect is replaced 
by two unlinked Saturn rings. Thereafter the colloids are 
displaced by competing forces of about w 30pN resulting 
from the topological restructuring of the disclination pat¬ 
terns following the field modulation, as detailed in SM. 
In particular, on one hand, the switching on of the field 
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FIG. 3: An equilibrated x-aligned cholesteric dimer is driven 
to the orthogonal, z-stacked configuration by a cyclic electric 
field modulated as a square-wave. 

favours the approach of the colloids and their alignment 
perpendicular to the field. On the other hand, whenever 
the field is switched off the colloids are pushed further 
apart and towards a vertical alignment. At steady state, 
the balance of these competing effects leads the colloids 
to stack approximately on top of each other thus forming 
an columnar arrangement along the z direction, as shown 
in Fig. [3] It should be emphasized that the observed col¬ 
loidal reorientation and ordering along the field direction 
is a purely non-equilibrium effect driven by the defect re¬ 
arrangement: accordingly, the free energy is higher for 
the columnar arrangement than for the figure-of-eight 
disclination. Importantly, the vertical colloid alignment 
is robustly attained upon varying the system size or by 
replacing the square-wave modulation with a sinusoidal 
one (see SM [36]). 

We finally turn to the case of several interacting par¬ 
ticles, which we address by considering a colloid suspen¬ 
sion of about 5% in volume fraction. Starting from an 
equilibrated random initial condition (Ai), we find that, 
as in the case of dimers, a time-dependent field (maxi¬ 
mal strength £ 3 , along the z direction) drives a collective 
alignment and clustering of the particles into colloidal 
columns along the field direction. This non-equilibrium 
ordering is essentially driven by the rearrangement of the 
disclination network created after field removal, that can 
displace particles efficiently due to impulsive forces simi¬ 
lar to those which we observed in the case of a dimer. At 
subsequent switchings on of the field, the defect-mediated 
interactions induce tighter columnar packing. This as¬ 
sembling and ordering process can be controlled by tun¬ 
ing not only the field strength or frequency, but also 
its direction. For instance, starting from the state of 
Fig. 0 (A 2 ) and then cycling the field along x and y, the 
colloidal columns rearrange into more disordered aggre¬ 
gates, which mainly lie on the xy plane (Fig. a (a*)). 

Importantly, the outcome is completely different when 
a static field of same strength and direction is applied 
(see Fig. [4] bottom panel). In this case, the kinetics en¬ 
tails a fast transition to the nematic state, where the 
colloids are accompanied by Saturn rings or figure-of- 
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FIG. 4: Top Panel: Controlled columnar and planar align¬ 
ment of several colloids with time-dependent electric fields. 
An equilibrated random dispersion of colloids (Ai) is first 
subjected to a cyclic field along z, resulting in vertical colum¬ 
nar alignments (A 2 ). Next, the repeated on-off switching of 
the field in the x, and y directions rearranges the colloids in 
horizontal planes (A 3 ). Bottom Panel: the same initial config¬ 
uration (Ai) is only slightly perturbed (Bi) by the application 
of a static electric field in the z axis. 


eight defects, but we observe negligible particle motion, 
so that the colloidal dispersion this time remains disor¬ 
dered. In conclusion, we have studied the physics of col¬ 
loidal dispersions in cholesterics under a time-dependent 
electric field. For a single particle, we showed that by 
tuning either the strength, or the time modulation of 
the field, it is possible to reach a variety of metastable 
states, characterised by distinctively different disclina¬ 
tion networks. Colloidal dimers and dispersions can also 
be repositioned under a time-dependent field. Most re¬ 
markably, a time-dependent field can rearrange an ini¬ 
tially random dispersion into a set of field-aligned col¬ 
loidal columns. This global repositioning is possible be¬ 
cause field-induced changes of the medium symmetry 
generate transient defects which, in turn, provide im¬ 
pulsive nonequilibrium forces favouring the axial align¬ 
ment of neighbouring colloids. We note that this ef¬ 
fect is robust upon replacing the dynamical evolution of 
eq. [2] with a more detailed one where hydrodynamic ef¬ 
fects [40j are accounted for with a hybrid Lattice Boltz¬ 
mann approach |41j . As detailed in SM, we carried out 
this computationally-intensive analysis for the dimer case 
and observed the same two previously-described mecha¬ 
nisms which compete for repositioning of the colloids. In 
particular: at field on, colloids tend to be drawn closer 
and towards an alignment in the xy plane while, at field 
off, they are pushed apart towards a 2 alignment. The 
cyclic variations of distance and vertical alignment, which 
are appreciably larger with hydrodynamics than without, 
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can be controlled by varying the interplay of these com¬ 
peting forces, e.g. with the relative duration of the field 
on and off intervals, see SM. We hope that such rich phe¬ 
nomenology can stimulate further theoretical and experi¬ 
mental studies on colloid-liquid crystal composites under 
time-dependent fields, or other nonequilibrium perturba¬ 
tions. 
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